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Abstract 



We study shear viscosity in weakly coupled hot pure gauge field QCD 
theory basing on transport theory and the Kubo formula using the closed 
time path formalism (CTP) of real time finite temperature field theory. We 
show that the viscosity can be obtained as the integral of a retarded three- 
point function. Non-perturbative corrections to the bare one loop result can 
be obtained by solving a Schwinger-Dyson type integral equation for this 
vertex. This integral equation represents the resummation of an infinite series 
of ladder diagrams which all contribute to the leading-log order result. We 
show that this integral equation has exactly the same form as the linearized 
Boltzmann equation and explain the reason behind this formal equality. 

PACS numbers: ll.lOWx, ll.lSTk, 11.55Fv 
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I. INTRODUCTION 



Fluctuations occur in a system slight perturbed away from equilibrium. The responses 
to these fluctuations are described by transport coefficients which characterize the dynamics 
of long wavelength, low frequency fluctuations in the medium [1,2]. The investigation of 
transport coefficients in high temperature gauge theories is important in cosmological ap- 
plications such as electroweak baryogenesis [3], and in the context of heavy ion collisions 
[4,5]. 

There are two basic methods to calculate transport coefficients: transport theory and 
response theory [6-12] . Using the transport theory method one starts from a local equilbrium 
form for the distribution function and performs an expansion in the gradient of the four 
velocity field. The coefficients of this expansion are determined from the classical Boltzman 
equation [12,13]. In the response theory approach one divides the Hamiltonian into a bare 
piece and a perturbative piece that is linear in the gradient of the velocity field. One 
uses standard perturbation theory to obtain the Kubo formula for the viscosity in terms of 
retarded Green functions [7,13]. These Green functions are then evaluated using equilbrium 
quantum field theory. As is typical in finite temperature field theory, it is not sufficient to 
calculate perturbatively in the coupling constant even for the weak coupling case. There are 
certain infinite sets of diagrams that contribute at the same order in perturbation theory 
and have to be resummed [13-15]. 

In this paper, we want to compare these two methods. The response theory approach 
allows us to calculate transport coefficients from first principles using the well understood 
methods of quantum field theory. On the other hand, the transport theory approach involves 
the use of the Boltzman equation which is itself derived from some more fundamental theory 
using, among other things, the quasiparticle approximation. In this sense, the response the- 
ory approach is more fundamental than the transport theory method. However, the response 
theory approach can be quite difficult to implement, even for a high temperature weakly 
coupled scalar theory, because of the need to resum infinite sets of diagrams [13,14]. These 
considerations motivate us to understand more precisely the connection between the more 
practical transport theory method, and the more fundamental response theory approach. 

Some progress has already been made in this direction. It has been shown that keeping 
only terms which are linear in the gradient expansion in the transport theory calculation is 
equivalent to using the linear response approximation to obtain the usual Kubo formula for 
the shear viscosity in terms of a retarded two-point function, and calculating that two-point 
function using standard equilibrium quantum field theory techniques to resum an infinite 
set of ladder diagrams for scalar fields [13,15-17]. This result is not surprising since it has 
been known for some time that ladder diagrams give large contributions to n-point functions 
with ultra soft external lines [18-22]. Some studies on transport coefficients and ladder re- 
summation for hot gauge theory using Kubo formule within the imaginary time formalism 
were reported in [23-25] . Due to the need for resummation, the diagrammatic approach to 
transport coefficients in hot field theories, based on Kubo formulae, has developed a rep- 
utation for being extremely cumbersome. This is particularly true for the ITF formalism 
with its need for analytical continuation and the corresponding complicated cutting rules 
[24] As a result, an alternative approach based on the Boltzmann equation in kinetic theory 
has recently become more popular [12] .With this appraoch, the transport cofficients in high 



2 



temperature theory has recently been studied up to leading order in [12]. But these calcula- 
tions are not trivial either, and they cannot fully replace the kind of intuitive insight into the 
underlying physical mechanisms which is provided by a diagrammatic analysis. Our goal is 
therefore to calculate the shear viscosity of hot gauge theory from both transport method 
and diagrammatic method in real time formalism , and thus understand the connection of 
these two methods. 

This paper will be organized as follows. In section II we define shear viscosity using 
a hydrodynamic expansion of the energy-momentum tensor. In section III we calculate 
shear viscosity using the transport theory method by performing a gradient expansion on 
the Boltzmann equation and obtaining a linearized Boltzmann equation. In section IV we 
calculate the shear viscosity using the Kubo formula which relates the shear viscosity to 
the two-point retarded Green function of the viscous-shear stress tensor. Starting from the 
Kubo formula we calculate the shear viscosity at leading-log approximation using standard 
techniques of finite temperature quantum field theory in real time formalism. We show that 
the shear viscosiy can be obtained as an integral over a three-point vertex. This three-point 
vertex satisfies an integral equation, which represents the resummation of an infinite series 
of ladder diagrams. We show that this integral equation representating ladder resummations 
has the same form as the linearized Boltzmann equation and explain the reason behind this 
formal equality. We discuss our results and present our conclusions section V. 



II. VISCOSITY 

In a system that is out of equilibrium, the existence of gradients in thermodynamic 
parameters like the temperature and the four dimensional velocity field give rise to ther- 
modynamic forces. These thermodynamic forces lead to deviations from the equilibrium 
expectation value of the energy momentum tensor which are characterized by transport 
coefficients like the thermal conductivity and the shear and bulk viscosities. In order to 
separate these different physical processes we decompose the energy-momentum tensor as, 

T^''' = eu^u" - pA^^ + P^'u" + P^u" + t^"" ; A^^ = g^^ - u^u^. (2.1) 

The quantities e, p, P^ and vr^j, have the physical meanings of internal energy density, 
pressure, heat current and viscous shear stress, respectively. The four vector u^{x) is the 
four dimensional velocity field which satisfies u^{x)u^{x) = 1. The expansion coefficients 
are given by 

P^ = A^^u^T^^, (2.2) 
TT^, = (A^,A,^ - ^A,^A,^)r"^. 

The viscosity is obtained from the expectation value of the viscous shear stress part of the 
energy momentum tensor. We expand in gradients of the velocity field and write. 
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= df^u, + d,u, - '^^^,^p^^Pu'' (2.3) 

where 77 is the coefficient of the term that is hnear in the gradient of the four velocity, defined 
as shear viscosity. 

Throughout this paper we work with weak-coupled SU{?>) Yang-Mills theory. The La- 
grangian for this theory is 

F;, = d,Al - d^Al - igr^'A^Al (2.4) 

where the coupling constant -C 1. In local rest frame, the traceless spatial part of the 
energy-momentum tensor is 



III. SHEAR VISCOSITY FROM KINETIC THEORY 

Kinetic theory and the Boltzman equation can be used to study transport properties 
of dilute many-body systems. One assumes that, except during brief collisions, the system 
can be considered as being composed of classical particles with well defined position, energy 
and momentum. This picture is valid when the mean free path is large compared with the 
Compton wavelength of the particles. At high temperature, the typical mean free path of 
thermal excitations in hot QCD is 0{l/g'^T\ng~^) and is always larger than the typical 
Compton wavelength of effective thermal oscillations which is 0{l/gT) . We introduce a 
phase space distribution function f{x,k) which describes the evolution of the phase space 
probability density for the fundamental particles comprising a fluid. In this expression and 
in the following equations the underlined momenta are on shell. 

The Boltzmann equation of pure gauge field QCD describing the evolution of the gluon 
distribution function f{x, k) has the form: 



d d ^ d 



/(k,x,i) = C[/]. (3.1) 



The external force Fext term will only be relevant in discussing the electrical conductivity 
. For calculations to leading logarithmic order in g, and for the transport coefficient un- 
der consideration, it will be sufficient to include in the collision term C[f] only two-body 
scattering processes, so that 

C[f] = l I C?ri2^3fc[/l/2(l + /3)(l+/fc)-(l + /l)(l + /2)/3A] (3.2) 
Z J 123 

with /j := f{x,p^), fk '■— fixjk). The symbol dri2^3k represents the differential transition 
rate for particles of momentum Pi and P2 to scatter into momenta P3 and K and is given 

by 
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where T is the multiparticle scattering amphtude. 
The form for f{x, k) in local equilibrium is, 

/^°^ - ,,i.)u,L^ _ 1 ; A^.:=l + 2n, (3.4) 

We study the Boltzmann equation in the hydrodynamic regime where we consider times 
which are long compared to the mean free time and describe the relaxation of the system 
in terms of long wavelength fluctuations in locally conserved quantities. For a simple fluid 
without any additional conserved charges, the only locally conserved quantities are energy 
and momentum. To solve the Boltzman equation in this near equilibrium hydrodynamic 
regime, we expand the distribution function around the local equilibrium form using a 
gradient expansion. We go to a local rest frame in which we can write u{x) — 0. Note that 
this does not imply that gradients of the form diUj must be zero. In the local rest frame 
(2.3) becomes, 

H'^ij) = -^Hij H (3.5) 

2 

Hij = diUj + djUi - -6ij{d ■ u) 

In all of the following expressions we keep only linear terms that contain one power of H^j 
and quadratic terms that contain two powers of H-ij, since these are the only terms that 
contribute to the shear viscosity we are trying to calculate. 
We write, 

/ = + f^'^ + • • • (3.6) 
with, /(I) Using (3.4) we obtain, 

A;^9^/(o) = /3nfc(l + nfc) \hAk)H^, 

/(^) -nk{l + nk)(t>k ; <t>k = P^Bij{k)Hij (3.7) 



where we define 



and write. 



iij{k) = ikkj - ^5ij) ; Iij{k) = eiij{k) (3.8) 



B{k)ij = iij{k)B{k) (3.9) 
The viscous shear stress part of the energy momentum tensor is given by 

r rl^k 1 



5 



Using the expansion (3.6) and rotation invariance 

kikjB{k) yijk^B{k) 

kikjkikmB{k) j^iSijSim + SuSjm + SimSji)k'^B{k) (3.11) 

we obtain the hnear contribution from (3.7) 

The lowest order term gives zero. Comparing with (3.5) we have, 

Thus we have shown that the shear viscosity can be obtained from the function B{k) . 
This function are the coefficient of the hnear term in the gradient expansion of the distri- 
bution function. Next we will show that this functions can be obtained from the linearized 
Boltzmann equation. 

Using the expansion (3.6) and (3.7) one can linearize Boltzmann equation 

ra^/°(x,^)=C[/(°);/W] (3.14) 

where we keep terms hnear in f^^^ on the right hand side. The left hand side gives, 

k!'di,f{x,k) = -nk{l + nk)P^Iij{k)Hij (3.15) 

The right hand side is 

C[f^'^]^-ll dru^3k{l + ni){l+n2)nsnk{(pk + (l>3-(pi-h) (3.16) 

Z J123 

Using the definition of given in (3.7) and comparing the coefficients of Hij on both side 
of the Boltzman equation we obtain, 

I,,{k) = \j dr,2^3fe ^^ + ^;^|^ + ^^^^^ [i?.,(p3) + B,,{k) - B,,{pJ - B,,{pJ ] (3.17) 

J 123 J- ~n '^k 

This result is an inhomogeneous linear integral equation which can be solved self-consistently 
to obtain the function Bij{k). 



IV. VISCOSITY AND LADDER RESUMMATION IN REAL TIME FORMALSIM 

The Kubo formulae allow us to use quantum field theory to calculate nonequilibrium 
transport coefficients. The results should be the same as those obtained in the previous 
section using transport theory. The Kubo formula expresses the shear viscosity in terms of 
a retarded two-point Green function [7,13,15] 
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v = TR:rr^^iH^R^^^]Uo=o (4.1) 

where Gr is the retarded two-point Green function of the viscous shear stress 

Gr{x, t; x', t") = -ie{t - f)[Kij{x, t),TTij{x', t")] (4.2) 

We obtain a diagramatic expansion for the viscosity from (4.1) . We use the closed time 
path formulation of finite temperature field theory, and work in the Keldysh representation. 
Several reviews of this technique are available in the literature [27-31]. The closed time path 
integration contour involves two branches, one running from minus infinity to positive infinity 
just above the real axis, and one running back from positive infinity to negative infinity just 
below the real axis [29,31]. All fields can take values on either branch of the contour and 
thus there is a doubling of the number of degrees of freedom. It is straightforward to show 
that this doubling of degrees of freedom is necessary to obtain finite Green functions. 

The two-point function or the propagator can be written as a 2 x 2 matrix of the form 

where Du is the propagator for fields moving along Ci, D12 is the propagator for fields 
moving from Ci to C2, etc. 

The retarded and advanced propagators are given by the combinations, 

Dr = L>n - D,2 

Da = Du - D21 . (4.4) 

The IPl part of the two-point function, or the polarization insertion, is obtained by trun- 
cating legs. The retarded and advanced parts are given by, 

Ur = Hn + 

= Hn + . (4.5) 

The situation is similar for higher n-point functions. For example, the three-point func- 
tion which is retarded with respect to the first leg is given by 

Tri = Till + Tn2 + Tm + Tu2 ■ (4.6) 

The other three-point vertices that we will need are: 

^R2 = Till + rii2 + r2ii + ^212 

Trs = Till + ri2i + r2ii + r22i 
— Till + ri2i + r2i2 + r222 

In our case, we want to obtain a perturbative expansion for the correlation functions 
of composite operators Dr{X,Y) which appear in (4.1). We define the vertices T^j as the 
vertice obtained by truncating external legs from the following connected vertice : 
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rgS? = {nn,,{X)Al{Y)AUZ)) (4.7) 

where is the operator that time orders along the closed time path contour. From this 
definition, one can draw the skeleton diagram of the two-point correlation function of the 
the viscous shear stress tensors with a full vertice and full propagators as 




Fig. 1: Skelon diagram of the two-point correlation function for shear viscosity from 
linear response. The dashed external line represents the composite operator TTij 

The coupling vertex between the operator tt*-' and two gluons with incoming momenta 
P, K and indices (/x, a), (z/, h) respectively can be read 

- (p'K^'S^'' + p'K'^S''' + P'^k'S^f' + P^k^'^) - '^d'^ {P-'k^'S'^' + K^'p'^S''")], (4.8) 

where P = {po, p), P ■ K — poko — p ■ k and po is the energy of the corresponding gluon. 
For g = 0,P = -K and on shell gluons = 0, r{K, Q, P) = r{K), we have: 

P^AK)T:^!:!2{K) = -2r^/.,(A:), 
PliK)r^SiK)PjpiK) = -2d^%{k)P^AK) (4.9) 

where Pi^^,{K) is the transverse projector defined as Pqq = P^^ = 0, and Pfjik) = 6ij — kikj. 
kj represents unit vector of k. Rotation invariance leads for the full vertice: 

PUK)r^SMK)Pjfs{K) = -25^'hj{k)P^^{K)r{K) (4.10) 

Since for the hard gluons, only the transverse part dominates the leading-log contribution, 
we can replace the full gluon propagtor by its transverse part. 

D''^'{K) ^ P^^DiK) (4.11) 

with 
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where IIt{K) represents the transverse gluon self-energy. 

Using the Feynman rules of the closed time path formalism (CTP) of real time finite tem- 
perature field theory,we write the two-point correlation function of the composite operators 



TT'i j as 



GM = i J dKVtTiK, g, -K - Q)iDZ{K)iD^,^{K + g)r(°)^^>^(p)r, (4.13) 

shown diagramatically in Fig.l. Where / dK = / , and the indices {a,b,c,d} are 
Keldysh indices and take values {1,2}. The gloun propagators here are full resummed 
propagators. We perform the sum over Keldysh indices using the Mathematica program 
described in [32]. 

After performing the sums over Keldysh indices and making use of Eqs (4. 9), (4. 10), we 
find imaginary part of the retarded two-point function 

ImGniQ) ^ -\ j dKk\Nk+g - N^^R^iK, Q, -K - Q)Da{K)Dr{K + Q) + h.c] (4.14) 

In obtaining this result we have used the fact that we will eventually take the limit Q ^ Q 
to obtain the viscosity (4.1). This hmit gives rise to what is known as the pinch effect: using 
an obvious notation we write the pairs of propagators Da{K)Dji{K -\- Q) as akrk+q, terms 
with a product of factors a^rk+q and r^afc+y contain an extra factor in the denominator that 
is proportional to the imaginary part of the sclf-cnergy relative to terms with products of 
propagators a^a^j^q or r^Tk+q- Thus terms proportional to a^ak^q or r^Tk+q can be dropped. 
This is the so-called pinch hmit, because the large terms occur when the contour is "pinched" 
between the poles of the two propagators, which gives rise to this enhencement factor in the 
denominator that is proportional to the imaginary part of the self-enegy. We regulate the 
pinching singularity with the imaginary part of the self-energy and obtain [15], 

rfcOfc+g - ^j^^j-^T ' Pk = K^k - ak) (4.15) 

where is the retarded part of the transverse gluon self-energy. 

Now we expand in go and keep only the term proportional to Qq since this term is the 
only one that contibutes to (4.1) . There is an over all thermal factor of the form Nk — Nk+q 
in (4.14). The expansion of this thermal factor is straightforward: 

Nk - Nk+q = 2go/3nfe(l + n^) + • • • (4.16) 

Consider the behaviour of the vertices when — > 0. We introduce the notation 
T{K,Q,—K) ■=T[K) . It is straightforward to show that 

Vr2{K, Q, -K-Q)^ RerR2{K) (4.17) 

Using these results to simplify (4.14) and substituting into (4.1) we obtain, 
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dKk'^ PkUkil + Uk) 



RerR2(K) 

imnj 



(4.18) 



Comparing with (3.13) we see that the results are identical if we use (3.9)and identify 



where momentum K must be on new mass shell: SlK"^ — Rell^) 

It has been known for some time that in scalar theory, the set of diagrams which give the 
dominate contributions to the viscosity are the ladder diagrams. These diagrams contribute 
to the viscosity to the same order in perturbation theory as the bare one loop graph, and 
thus need to be resummed [13,26,15]. Following the similar power counting rule one can 
show that the set of ladder diagrams shown in Fig. 2 give the dominate contributions to 
the viscosity in pure gauge field QCD theory at high temperature. This effect occurs for 
the following reason. It appears that the ladder graphs are suppressed relative to the bare 
vertex by extra powers of the coupling, which come from the extra vertex factors that one 
obtains when one adds rungs (vertical lines). However, these extra factors of the coupling 
are compensated for by a kinematical factor. This factor arises through the pinch effect 
which is described above. The addition of an additional rung in a ladder graph always 
produces an extra pair of propagators of the form akVk+q or rkdk+q- Products of this form 
contribute a factor which produces an enhancement. This factor occurs when the contour 
is "pinched" between the poles of the two propagators which gives rise to a contribution in 
the denominator that is proportional to the imaginary part of the self-energy. 

In order to include those ladder diagrams, which dominate the contribution to the vis- 
cosity, we obtain the vertex Tij{P, Q, —P — Q) as the solution to the integral equation shown 
in Fig.2. 



Fig.2: Integral equation for ladder resummation, the blob stands for the full vertex 
r . Thick lines represent for hard transverse gluons, where a finite width is needed 
to regularize the pinch sigularity. While thin lines represent soft gluons, where HTLs 
resummed propagators are required. 

Since we only calculate up to the leading-log order, in order to simplify the calculations, 
several points are in order: 

1) We consider weak coupled degrees of freedom {g <ti 1); 

2) Hard degrees of freedom with the momentum of order T for the rails dominate the 
contribution; 




(4.19) 




K 



K+Q 
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3)Small momentum transfer {S ^ gT),t chanel dominates the contribution. Namely, the 
momenta on the nmgs (vertical lines) are soft of order gT. 

In order to to get rid of color indices and Lorentz indices, We contract 5"'^P^j^{K) with 
both sides of above equation and make use of 

^ab ^aa' ^bb' j^a'cd ^ce jefg ^ff ^gg' jf'e-'g' ^e'c' jb'd'c' 

— jaed jefg jaed jefg _ 20) 

and 

PUR) 7a,«,a' Pl5> (R) 15,.', 5' = (4.21) 

where ia,f3,j{P: T, K) — (P—T)^ 6ai3+{T—K)a Sf)j+{K—P)pSa-y comes from the triple-gluon 
vertex. 

After some algbra calculations, we obtain a matrix integral equation in Keldysh space 

r[Z{K, Q, -K-Q)^ r(°)''- - UNl j dP dR T'ZiP, Q, -P - Q)D^{P)D,c{P + Q) 
■D^{R + K- P)T^D,a{R)PrR.D-''{T)P,'Rn'D^'^'{T)T^ (4.22) 

We perform the sums over the Keldysh indices using the Mathematica program in [32] 
and simplify the result by taking Q to zero, keeping only the pinching terms, and using 
(4.15). We obtain, 

r^2(^) = kmk - -S^ik^ + ^ / dPdRdR! {2%f5\P + R-R'-K) 

|r|VfiPfi' j^gT(p) (l + + nR)nR>/{l + Uk) (4.23) 

where = IGg'^lPrrg'^Rnl'^ , which is recognized as the matrix element squared for the 
collision of the hard gluon in t chanel as shown in Fig. 3(a). Here rj'* is the hard thermal 
loops resummed retarded propagator for the soft transfer momentum S. 




Fig. 3: (a) Two-body giuon-ghion clastic scattering in t chanel. (b) Sclf-cnergy 
describing collisions in the (resummed) Born approxximation. The thick lines refer 
to the the hard colliding gluons. The thin lines with a blob denote soft gluon 
propagators dressed by the HTL 
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Note that this integral equation is decoupled: the only three point vertex that appears 
is r^2- To simplify this expression further we write Tij = Iij{p)r{P), use (4.19) and sym- 
metrize the integral on the right hand side over the three integration variables P, R' and R. 
We multiply and divide the left hand side by Imn^(P) and replace this expression in the 
numerator by the HTL result [21,33,34] 

Imn^(X) = \ (^^) / dPdR dR' {2t:Y8\P + R-R'- K)\T\'' 

PpPrPR'i^ + np){l + nR')nR (4.24) 
Rearranging we obtain [15], 

I^k,k)i^ = — [ dP dR dR' {2n)%P + R-R'- K)\T\^ + ""f + ^"^'^""^ (4.25) 
2 J (1 + Uk) 

PpPrPR'[Blm{R') + Bim{K) - BUP) - BUR)] 

where we have used (4.17) and (4.19). When the delta functions are used to do the 
frequency integrals, this equation has exactly the same form as the equation obtained from 
the linearized effective Boltzman equation(3.17), with a new mass shell describing thermal 
excitations. Comparing (3.13) and (3.17) with (4.18) and (4.25) we conclude that calculating 
shear viscosity of the hot pure gauge field QCD using effective transport theory by keeping 
only first order terms in the gradient expansion, is equivalent to using the Kubo formula 
obtained from response theory, with a three point vertex obtained by resumming ladder 
graphs. 

Now we have established an interesting connection between kinetic theory and diagram- 
matic linear response theory for the nonabilian gauge theory. We obtain formally the same 
expressions for the viscosity in both approaches. It is desirable to understand the deeper 
reasons behind this formal equality. 

The viscosity is obtained from the expectation value of the viscous shear stress tensor. 
There are two approaches to calculate the expectation value. The first method is from 
kinetic theory. With this approach, we use a distribution function / to evaluate the the 
expectation value of the viscous shear stress tensor(see Eq.(3.10)). This distribution function 
is expanded around a local equlibrium function f^^'^ up to the first order /*^^) of the gradient 
expansion . The second approach is by linear response theory. The expectation value 
is calculated from the non-equilibrium density matrix p [7, 17], which can be expanded 
around the local equilibrium density matrix po up to the linear response order in gradient 
expansion , p = po + The expectation value of a physics quantity does not depend on the 
approaches one uses, either from the kinetic theory by the distribution function , or from 
the linear response theory by the density matrix under the same approximation in gradient 
expansion. Therefore B{k) from (3.13) and from (4.19)should be the same although they 
seem to describe completely two different things! The former one comes from the first- 
order distribution function Z*^^); The later, though related to a retarded 3-point function, is 
resulted from the first-order matrix density pi. Physicaly , they both describe how far the 
system is perturbated away from the local equilibrium. In fact, both methods describe the 
first-order deviation of the system away from the local equilibrium, one of them using the 
description of the distribution function, while the other one using the density matrix. 
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On the other hand, both appraochcs also contain the same physics processes (dia- 
gramms), namely two body collisions. In order to understand this point, let's recall that the 
Boltzmann equation can be derived from the Schwinger-Dyson equation of the full 2-point 
functions in real-time formalism under certain approximations , like Wigner transformation 
and truncation schemes etc [21]. Different truncation schemes (or self-energy insertions ) 
lead to different collision terms. The Boltzmann equation we used only contains 2 — >^ 2 colli- 
sion terms, which come from the cut of the self-energy diagrams containing 2^2 processes 
as shown in Fig. 3 . The cutting rules tell us that the cuts of self-energy diagramms give 
physics scattering amplitudes. The diagrammatic method from the integral equation of the 
full vertex is essentially like a resummation for the full 2-point correlation diagram of the 
stress tensor. Again in this resummation, only the diagramms containing 2 — > 2 processes 
are included. 

Because of these reasons above, we expect these two approaches should be able to give 
the same result for any theories under the same approximation. 

V. SUMMARY AND CONCLUSIONS 

We have studied shear viscosity of the nonabelian pure gauge theory using two different 
methods. The first method uses standard transport theory. We start from a local equilibrium 
form for the distribution function and perform a gradient expansion. We calculate the shear 
viscosity by expanding the Boltzman equation and obtaining a linearized Bolzman equation 
that can be solved consistently. The second technique uses the Kubo formula for linear 
response, which allows us to calculate the shear viscosity from the retarded two-point Green 
function of the viscous shear stress tensor. The transport theory calculation involves the use 
of the Boltzman equation, which is itself obtained from some more fundamental theory. The 
response theory calculation uses the well known methods of equilibrium finite temperature 
quantum field theory and is, in this sense, more fundamental. However, the calculation is 
complicated by the need to resum infinite sets of diagrams at finite temperature. At leading 
order, it is well known that a correct calculation of the linear response coefficient in scalar 
theory involves the resummation of ladder graphs. 

We have identified precisely which diagrams in hot pure gauge QCD need to be resummed 
by studying the connection between the transport theory calculation and the response theory 
calculation in real time formalism. We have shown that calculating shear viscosity using 
transport theory by keeping terms that are linear in the gradient of the velocity field in 
the expansion of the Bolzman equation, is equivalent to calculating the same shear viscosity 
from quantum field theory at finite temperature using the linear response Kubo formula with 
a vertex given by a specific integral equation. This integral equation is the same as that 
recently found by Arnold, Moore and Yaffe [12] by analyzing the infrared divergences of the 
linearized collision integrals without screened interactions. This integral equation is also the 
same as the linearized Bolzmann equation derived by Blaizot and lancu from the Kadanof- 
Baym equations by using a gradient expansion and keeping only linear terms [21]. Thus, this 
fact constitutes a non trivial check of the formalism we have used. This integral equation 
shows that the complete set of diagrams that need to be resummed includes infinite ladder 
graphs. We explain the reasons behind this formal equality between this two approaches. 
First, both approaches describe the first-order deviation of the system away from the local 
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equilibrium, one of them using the description of the distribution function, while the other 
one using the density matrix. Secondly both methods involve the same physics processes or 
diagramms ( two body collisions). Therefore, We expect this interesting connection between 
these two appraoches should exist for any other theories under the same approximations. 

There are several directions for future work. The study here is only valid up to leading- 
log approximation of the transport coefficients. Going beyond to this approximation, one 
needs to include certain inelastic LPM suppressed splitting processes and thus needs to 
resum more complicated diagrams [12,35]. It would be interesting to generalize this work to 
the full leading order of of the shear viscosity with real QCD with quark flavors in real time 
formalism . 
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